Memory function formalism approach to electrical conductivity and optical response of dilute 

magnetic semiconductors 
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A combination of the memory function formalism and time-dependent density-functional theory is applied 
to transport in dilute magnetic semiconductors. The approach considers spin and charge disorder and electron- 
electron interaction on an equal footing. Within the weak disorder limit and using a simple parabolic approxi- 
mation for the valence band we show that Coulomb and exchange scattering contributions to the resistivity in 
Gai_a;Mna; As are of the same order of magnitude. The positional correlations of defects result in a significant 
increase of Coulomb scattering, while the suppression of localized spin fluctuations in the ferromagnetic phase 
contributes substantially to the experimentally observed drop of resistivity below Tc- A proper treatment of 
dynamical screening and collective excitations is essential for an accurate description of infrared absorption. 
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The idea of utilizing the caiTier spins in new electronic de- 
vices provides the basis for the rapidly developing field of 
spintronics [7]. A unique combination of magnetic and semi- 
conducting properties makes dilute magnetic semiconductors 
(DMSs) attractive for various spintronics appUcations 10]. A 
lot of attention is drawn to Gai_3;Mnj;As since the discovery 
of its relatively high ferromagnetic transition temperature |0], 
with a current record of Tc = 159 K^. 

The sensitivity of magnetic and transport properties of 
Gai_2:Mna;As to the details of growth conditions |4] and post- 
growth annealing ifslliil?] points to the crucial role played by 
the defects and their configuration, and has stimulated intense 
research on the structure of defects and their influence on the 
various properties of the system [8J. Most theoretical calcu- 
lations of the electrical conductivity and optical response in 
these materials, however, treat disorder within the relaxation 
time approximation, where the relaxation time is often treated 
as an adjustable phenomenological parameter. It is essen- 
tial, therefore, to develop a theory of electrical conductivity 
in DMSs with emphasis given to disorder 

Our approach is based on a combination of the memory 
function formalism 1^ [l^ [Till and time-dependent density- 
functional theory (TDDFT) |[l2|] . which not only goes beyond 
the simple relaxation time approximation for disorder scatter- 
ing, but allows one to consider key features of DMSs such as 
spin and charge disorder and electron-electron interaction on 
an equal footing. The present paper is thus an extension of our 
previous work 1 13], in which we first introduced the memory 
function formalism for transport in DMS, but did not address 
spin disorder and electronic many-body effects beyond static 
screening. 

Disorder in our model consists of the Coulomb potential 
of the charge defects and fluctuations of localized spins (the 
mean-field part of p-d exchange interaction enters the clean 
system Hamiltonian). Introducing a four-component nota- 
tion, one can express the disorder Hamiltonian as Hj = 

X]k^(^) ■ where the four-component disorder po- 



tential 
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is coupled to the four-component charge and spin density op- 
erator p''(k) = ■^EqErr'(c^'')rr' a+_k,r «q,r' • Here, ct'' 

(H — 1, +, — , z) is defined via the Pauli matrices, where cr^ is 
the 2 X 2 unit matrix, a"^ — {a^ ± iay )/2, and the sum in ([T]) 
runs over all defects. 

To describe transport in DMSs we employ the memory 
function formalism ^ [l^, [HI]. The central point of this ap- 
proach is the calculation of the current relaxation kernel (or 
memory function), whose imaginary part can be associated 
with the Drude relaxation rate. Our derivation of the memory 
function in spin- and charge-disordered media is based on an 
equation of motion approach for the current-cuiTent response 
function lfl4l [Tsll . Technical details of the derivation will be 
published elsewhere. Here we only give the final expression 
for the memory function in the long-wavelength limit in para- 
magnetic state: 



where 



XpMp^(k,t^) 



mnuj 



X 
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a,P — x,y,z (idi , n is the carrier concentration, averag- 
ing is performed over the magnetic subsystem Hamiltonian 
Hjn, and Xp'^p'^ (k, t^) are charge- and spin-density response 
functions associated with the operators p^. Eq. (|2]i contains 
the set of charge- and spin-density response functions of the 
disordered system and, strictly speaking, should be evaluated 
using a self-consistent procedure. This approach was real- 
ized in Ref. [17] to study a spin-independent system close to 
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FIG. 1 : Charge and spin relaxation rates for random and correlated 
impurities in Gao.95Mno.05As. See text for details. 



the metal-insulator transition. In our case, however, we as- 
sume that the disorder is weak enough so we can approximate 
Xpi^p" (k, oj) by their clean (disorder-free) system counterparts 
XpMp. (k,cj). 

The form of the disorder potential ([T]) allows us to sepa- 
rate in Eq. ^ contributions from charge and spin disorder 



M{uj) 



T„(a;) 



with 



= ^ kakp 
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^i(k) 



L(k,w), 



(4) 



and 



(5) 

The imaginary parts of Eqs. (|4| and (|5]l represent the charge 
and spin relaxation rates, which depend on frequency and, 
in general, on momentum. Calculations of the charge and 
spin relaxation rates have been performed for the case of 
Gao.95Mno.05As, assuming for simplicity parabolic disper- 
sion for holes. The material parameters used are: heavy hole 
effective mass m = 0.5 mo, dielectric constant e = 13, and 
exchange constant VJ = 55 meV nm'^, which corresponds to 
the widely used DMS p-d exchange constant Nof3 = 1.2 eV 

Most theoretical models for transport in DMSs assume ran- 
dom defect distributions. Monte-Carlo simulations of Timm 
et al. f]^, however, have shown that, driven by the Coulomb 
attraction, donor and acceptor defects in Gai_2;Mnj:As tend 
to form clusters. The main effect of such a clustering is 
ionic screening of the disorder Coulomb potential, which has 
been shown to be necessary to correctly reproduce the band 
gap, metal-insulator transition and shape of the magnetization 
curve. But the correlation of defect positions also gives rise 



to a momentum dependent impurity structure factor We here 
account for both effects. 

We consider systems with two types of defects, manganese 
ions in substitutional (Mnca) and interstitial (Mnj) positions. 
We treat the former as the acceptor centers that carry localized 
spins, and the latter as spinless double donors. The donor- 
acceptor cross term in Eq. ^ accounts for ionic screening. 

Correlation in the defects positions gives rise to a set of 
impurity structure factors in the product of the charge compo- 
nents of the disorder potential in Eq. (|4|: 



Ul{k) U{ (-k) = umUl (-k)-^ Su'{k). 



(6) 



Here, the structure factors can be expressed through the cor- 
responding pair distribution functions Pui , 
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^ / Pu' (R) cos(k • R)dR, (7) 
"0 Jv 



where i, i' label acceptors and donors and fio is the elementary 
cell volume. We use a simple model expression for Pai (for 
details, see fl3ll ). with parameters consistent with the Monte 
Carlo simulations of Ref. 

Positional correlations alone, however, are not sufficient to 
affect the spin relaxation rate (|5]l, orientational correlations of 
the spin fluctuations are also necessary. To account for these, 
we considered interacting spins described by the Heisenberg 

Hamiltonian H„i ~ '^j^j' "^n' i ' i' ^ ^^S^ 

temperature expansion to obtain the following expression for 
the product of the spin components of the disorder potential 
in Eq. (01: 



W^(k)Z^,(-k) 
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where ^x^v = x,y,z and the spin structure factor iSs(k) (ad- 
justed for orientational correlations) is 



Ss{k) = l- 



^SuniSMn + 1) 
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j>j' 



kuT 



cos (k • TLjj') 



(9) 

Here, Syin = 5/2 and Ni are the spin and number of localized 
moments, and the d — d exchange constants Jjj' are chosen 
to reproduce the typical experimental value of Tc = 150 K 
within the standard mean field approach. 

In Fig. [U we plot the static {lo = 0) relaxation rates (|4]i- 
(|5]l calculated for random (dashed lines) and correlated (solid 
lines) impurities in Gao.95Mno.05 As as a function of the level 
of compensation. Correlated impurities were modelled by 
clusters containing 10 Mnca with an average concentration 
Xc ~ 0.1 of substitutional Mn ions within the cluster. It is 
seen that the charge and spin contributions to the relaxation 
rate are of the similar magnitude. The combined effect of 
ionic screening and impurity structure factor results in a net 
increase of the charge relaxation rate for correlated impurities 
for the whole range of compensations. The increase is signifi- 
cant (up to 100%) and is sensitive to the cluster configurations. 
The latter might be controlled by the post growth annealing. 
The positional correlation of the scattering centers also leads 
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FIG. 2: Temperature dependence of the charge and spin relaxation 
rates for Gao.gsMno. 05 As with Tc — 150 K. Inset: charge relaxation 
rate, normalized to its value in the paramagnetic state. 



FIG. 3: Frequency dependence of the charge and spin relaxation rates 
for Gao.95Mno.05 As with electron-electron interaction taken into ac- 
count within static and dynamic RPA. 



to an increase of the spin relaxation rate for interacting spins 
at room temperature. This effect, however, is smaller than that 
for charge relaxation. 

Magnetic ordering is known to have a strong effect on the 
transport properties of DMSs [20]. The resistivity of opti- 
mally annealed samples reveals a pronounced maximum at 
critical temperature and decreases significantly for tempera- 
tures below Tc ^27]. Lopez-Sancho and Brey l'22'l proposed 
to explain the resistivity change in terms of the variation of 
the Fermi surface and the transport scattering time when go- 
ing from the paramagnetic to the feiTomagnetic phase. Their 
model, however, completely neglects the scattering off the 
fluctuations of localized spins. On the other hand, spin fluctu- 
ations are effectively suppressed in the ferromagnetic state. In 
the fully spin-polarized state scattering takes place only due 
to the quantum fluctuations of localized spins. 

In Fig. |2] we present the temperature dependence of spin 
and charge relaxation rates calculated according to Eqs. (|4|i 
and (|5]l for Gao.95Mno.05As with Tc = 150 K. Suppression 
of localized spins fluctuations below Tc leads to a significant 
reduction of the spin relaxation rate. Given the comparable 
magnitudes of charge and spin relaxations, the 70% drop in 
the latter translates into about 20% reduction in total resis- 
tivity, which is consistent with both experimental observa- 
tions and the calculations of Ref. |22]. This shows that the 
spin scattering is clearly not negligible in Gai_2;Mna;As, es- 
pecially if one considers effects associated with magnetic or- 
dering. 

We also determined the variation of the charge scattering 
rate with temperature, which was discussed in Ref. |22]. How- 
ever, we find its magnitude to be much smaller, see the inset 
in Fig. (|2]l. This is most likely due to the fact that in our cal- 
culations we used a simple parabolic band and isotropic spin 
splitting. Nevertheless, it is clear that further work is required 
to adequately describe both scattering mechanisms and their 
contributions to the experimentally observed drop in resistiv- 



ity. 

Most previous studies of (magneto)transport in DMS in- 
cluded electronic many-body effects only in the form of static 
Coulomb screening ||23[ I24ll25ll . However, this simplification 
ignores the role of dynamical many-body effects such as the 
coupling to plasmon modes. A major advantage of the mem- 
ory function formalism is that it allows to consider both dis- 
order and electron-electron interaction on equal footing. All 
carrier many-body effects in Eq. ^ including screening, cor- 
relations and collective excitations are absorbed in the set of 
density and spin-density response functions and, in principle, 
can be accounted for exactly by means of TDDFT [12]. In the 
original work by Gross and Kohn ll26ll . the interacting density- 
density response function of a homogeneous system is shown 
to be representable as 

4:7T 

X"^(q,w) = X(7\q,w) ^ /xc(q,^) , (10) 

with the exchange-correlation kernels f^c- 

Coupling to the charge plasmon mode already occurs on the 
level of dynamic RPA, corresponding to the first two terms on 
right hand side of Eq. ( fTOl l generalized for the multicomponent 
spin-charge density response functions. Fig. ^ shows the 
frequency dependence of charge and spin relaxation rates cal- 
culated for Gao.95Mno.05As within static and dynamic RPA. 
Coupling to plasmon modes results in a strong enhancement 
of the charge relaxation rate since it provides an efficient chan- 
nel to absorb the momentum from impurity scattering. 

This approximation, however, does not affect the spin re- 
laxation. To capture collective spin modes one has to go be- 
yond RPA and include exchange and correlation contributions 
in Eq. ( fTOl l. which can be done in the adiabatic local-density 
approximation ll27ll . However, our single band model does 
not produce any spin collective mode. In a a more realistic 
model with multiple valence bands, inter-valence band spin 
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collective modes may play a role. This is currently work in 

progress. 

To summarize, we have derived a general framework for 
combining the memory function formahsm with TDDFT in 
spin and charge disordered media, to study transport proper- 
ties of DMSs. Within the weak disorder limit and using a 
simple paraboUc model for the valence band we have shown 
that Coulomb and exchange scattering contributions to the re- 
sistivity in Gai_j;Mn2;As are of the same order of magnitude 
and should be taken into account simultaneously. The com- 
bined effect of ionic screening and impurity structure factor 
results in a net increase of the relaxation rate in systems with 
positional correlation of the defects. The suppression of local- 
ized spin fluctuations in ferromagnetic phase contributes sub- 
stantially to the drop of resistivity experimentally observed 
below Tc. Our calculations further suggest that the effects of 



collective electron dynamics on transport and optical conduc- 
tivity in DMS is significant. One can expect that there will be 
a distinct collective signature in the mid-infrared free-carrier 
absorption, which has not been included in most previous the- 
oretical studies. These results, in particular when validated in 
currently ongoing work which includes band structure details, 
should give valuable insight into the transport properties of 
DMSs. 
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